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Abstract 

In this paper, an exact smooth solution for the equations modeling the bedload transport of sediment in Shallow 
Water is presented. This solution is valid for a large family of sedimentation laws which are widely used in erosion 
modeling such as the Grass model or those of Meyer-Peter & Miiller. One of the main interest of this solution is 
the derivation of numerical benchmarks to valid the approximation methods. To cite this article: A. Namel, A. 
Name2, C. R. Acad. Set. Paris, Ser. I 340 (2005). 

Resume 

Une solution analytique du systeme de Saint- Venant couple a l'equation d'Exner. Ce papier presente 
une solution analytique pour le systeme modelisant le transport de sediments par le charriage. Cette solution est 
valable pour une grande famille de lois sedimentaires comme le modele de Grass ainsi que celui de Meyer-Peter 
& Miiller. Ce resultat est utile pour la validation des schemas numeriques. Pour citer eel article : A. Namel, A. 
Name2, C. R. Acad. Sci. Pans, Ser. I 340 (2005). 



1. Introduction 



Soil erosion is a consequence of the movements of sediments due to mechanical actions of flows. In the 
context of bedload transport, a mass conservation law, also called Exner equation [1], is used to update 
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the bed elevation. This equation is coupled with the shallow water equations describing the overland flows 
(see [2] and references therein) as follows: 

d t h + d x {hu) =0, (1) 
d t {hu) + d x (hu 2 + gh 2 /2) + ghd x z b = 0, (2) 
d t z b + d x q b = 0, (3) 

where h is the water depth, u the flow velocity, z b the thickness of sediment layer which can be modified 
by the fluid and g the acceleration due to gravity. The variable hu is also called water discharge and noted 
by q. Finally, q b is the volumetric bedload sediment transport rate. Its expressions are usually proposed 
for granular non-cohesive sediments and quantified empirically [3,4,5]. 

Many numerical schemes have been developed to solve system (1-3) (see [5] and references therein). The 
validation of such schemes by an analytical solution is a simple way to ensure their working. Nevertheless, 
analytical solutions are not proposed in the literature. Up to our knowledge, asymptotic solutions, derived 
by Hudson in [6], are in general adopted to perform some comparisons with approximated solutions. The 
solutions are derived for Grass model [3], i.e q b = A g v?, when the interaction constant A g is smaller than 
1CP 2 . In this paper, we propose a non obvious analytical solution in the steady state condition of flow. 



2. Solution of the equations 



In order to obtain an analytic solution, we consider q b as a function of the dimensionless bottom shear 
stress t b (see [5]). By using the friction law of Darcy-Weisbach, is given by 



T b = 



fu 2 



8(s-l)gd,' 



where / is the friction coefficient, s — p s /p the relative density of sediment in water and d s the diameter 
of sediment. The formulae of q b is usually expressed under the form 



(4) 



where r* r is the threshold for erosion, k an empirical coefficient and p an exponent which is often fixed 
to 3/2 in many applications. The expression (4) can be written in the simple form 



q b = Au 2 /, 

where the effective velocity u e and the interaction coefficient A are defined by 



(5) 



U = ^ T 



1 



l)gd 
1 p 



(6) 



8(a - l)gd 



Remark. The Grass model [3] is one of the simplest case by using p = 3/2, t* t = and an empirical 
coefficient A g instead of A. The Meyer-Peter & Miller model [4] is one of the most applied by using 
p = 3/2, k = 8, r* r = 0.047. The following result is valid for all models rewriting in form (5-6). 
Proposition 2.1 Assume that q b is defined by (4). For a given uniform discharge q such that t£ > r* r , 
system ( 1-3) has the following analytical unsteady solution 



2 



ax + (3 



A 



i/p 

1 

, h = q/u, 



_ u + 2gq 
^ z 6 = -at + z°. 



(7) 



where a, /3, C are constants and A, u cr are defined by (6). 

Proof. We are here concerned by the smooth solution. In view of the assumption hu = q 
(1-3) reduce to 



est, equations 



d t h = 0, 

d x (q 2 /h)+ghd x H = 0, (8) 
d t H + d x q b = 0, (9) 

where H = h + z is the free surface elevation. Differentiating equation (8) with respect to t and then 
equation (9) with respect to x, we obtain 

d xt H 



0. 



d 2 x q b = 0. 



(10) 



Note that we can write q b = qb(h,q) to have dtqb = dhqbdth + d q qbdtq = 0, so q b is not time-depending. 
Thank to (10), the expression of qb is obtained under the form 



q b = ax + f3, 

where a and (3 are constant. From (3), we obtain d t Zb = —d x qb 

z b = -at + z°{x). 
Moreover, from (5) we deduce the effective velocity as follows: 



-a to write 



(11) 
(12) 



ax + [3 



i/p 



Plugging (12) into the momentum equation (8) and using a direct calculation, we have 



e> x z h 



d x u 



z b 



2gq 



2gu 



+ C 



which concludes the proof. 

Remark. As h and u are stationary, the initial condition of (7) is (h, u, z b ). Moreover, the solution (h, u) 
applied to the Grass model is also an analytical solution of the Shallow Water Equations with the variable 
topography z°. Concerning the shallow- water model, other solutions can be found in [7]. 



3. Numerical experiments 

In this section, we consider the analytical solution (7) applied to the Grass model with q = 1, A g = 
a = (3 = 0.005 and C = 1. A relaxation solver is applied to approximate the solution of the model. We 
will not give here the details of the relaxation solver (for details see [8]), but just the relaxation model 
for the equations (1-3). Thus, we solve the following relaxation system: 



3 



J=500 



J=500 




x [m] x [m] 

Figure 1. Comparison between the exact solution and the relaxation method for : the water height and the topography (left) 
and the velocity (right). 



' d t h + d x (hu) = o, 

d t {hu) + d x (hu 2 + ?r) + ghd x z b = 0, 
a 2 

< d t n + u8 x tt + —d x u = 0, 
h 

d t Zb + d x q r = 0, 

d t q r + (jp - u 2 ^j d x z b + 2ud x q r = 0, 

that is completed with 7r = gh 2 jl and q r = q b at the equilibrium. Figure 1 presents the numerical result 
with J = 500 space cells, a CFL fix condition of 1 and T — 7s. We only notice little difference on the 
velocity, near the inflow boundary. 
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